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GRADIENT DYNAMICAL SYSTEMS ON OPEN SURFACES 
AND CRITICAL POINTS OF GREEN'S FUNCTIONS 

ALBERTO ENCISO AND DANIEL PERALTA-SALAS 



f^^ Abstract. We study the dynamics of the vector field on an open surface 

fvQ given by the gradient of a Green's function. This dynamical approach enables 

, us to show that this field induces an invariant decomposition of the surface as 

f^ the union of a disk and a 1-skeleton that encodes the topology of the surface. 

<We analyze the structure of this 1-skeleton, thereby obtaining, in particular, 
a topological upper bound for the number of critical points a Green's function 
|y^ can have. Connections between the dynamical properties of the gradient field 

and the conformal structure of the surface are also discussed. 



1. Introduction 



■3 

Cd Let M be a noncompact surface without boundary of class C°°, endowed with 

H a smooth complete Riemannian metric g. We denote by 

' ' g : (M X M)\ diag(M x M ) -> E 

►^ a Green's function of M, which is defined as a symmetric function (i.e., Q{x,y) = 

rvq G{y,x)) that satisfies the equation 

T-H for each y E M . That is to say, if one considers the action of the Laplace-Beltrami 

^^ operator of the manifold, Ag, on the Green's function g{x,y) (with respect to the 

f^ first variable x), one obtains a Dirac measure supported at the point y. 

" ' Our goal in this paper is to analyze the dynamical properties of the gradient of 

. . the Green's function. For this, we will find it notationally convenient to fix a point 

y G M, once and for all, and use the notation G := G{-,y) for the Green's function 
with pole y, which is smooth and harmonic in M\{y}. Therefore, the gradient field 
;J] we will study in this paper will be ^ gG. 

The study of Green's functions is a central topic in Riemannian geometry and 
geometric analysis. Hence, there is a vast related literature covering, among many 
other aspects, the existence of positive Green's functions [3l[TT], upper and lower 
bounds, gradient estimates and asymptotics [131 13], and the connection between 
Green's functions and the heat kernel [HI [TO] . 

What is somewhat surprising is that the dynamical properties of the gradient 
field VgG remain virtually unexplored, with the exception of the classical work 
of Brelot and Choquet |lj. Of course, a key issue in the study of this vector 
field is the analysis of the critical set of the Green's function. This question is 
of considerable interest by itself, and deeply related with other problems, set in 
significantly easier contexts, that date back to the 1950s (see e.g. [3D1 [T71 [TH [TO] 
and references therein). Indeed, some of these articles were motivated in part by 
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the fact that in Euchdean space the Green's function arises as the electric potential 
of a charged particle, so that its critical points correspond to equilibria and the 
trajectories of its gradient field are the force lines studied by Faraday and Maxwell 
in the XIX century (nontrivial contributions to this problem were made in the 
recent paper [7]). As a side remark, let us point out that another elliptic PDE 



(very different from ( 1.1 1) in which the analysis of the critical points of its solutions 
has recently attracted considerable attention can be found in [14^ . 

One reason why the study of the dynamical properties of the gradient field V^G 
(and in particular of the critical set of G) is hard to tackle, from the point of 
view of geometric analysis, is that the estimates for Green's functions obtained 
through PDE methods are not sufficiently fine to elucidate whether the gradient 
of G vanishes in a certain region. Moreover, the noncompactness of the underlying 
surface introduces complications related to the behavior of the Green's function at 
infinity. 

In this paper we will show how these difficulties can be overcome by exploiting 
the conformal properties of the surface and resorting to techniques of gradient 
dynamical systems. Our approach will lead to a topological upper bound for the 
number of critical points of the Green's function and a complete description of the 
local and global dynamics of the gradient field V^G. 

To some extent, the core of this paper is the remarkable heuristic principle we will 
now state, which links the dynamics of the gradient field ^ gG, defined in terms 
of solutions to an elliptic PDE, with the topology of the underlying surface. It 
must be stressed that this principle will be promoted to a rigorous statement (after 



introducing the necessary tools and notation) in Corollary 4.6 and Theorem 5.4 



Heuristic principle. Suppose that the surface M is of finite type. Then M can 
be decomposed as the union of a disk D and a (possibly disconnected) noncompact 
graph T , both of which are invariant under the local flow of the gradient field 'VgG. 
The disk consists of the pole y and the points of M whose ui-limit is y. The graph J- 
consists of the critical points of G, their stable components, and certain trajectories 
of VgG that escape to infinity. When suitably compactified, T is a connected graph 
that encodes the topology of the surface, the rank of the first homology group of J- 
being twice the genus of M . 

The characterization of the set J- that will emerge from the rigorous version of 
this heuristic principle yields, as a nontrivial application, the following topological 
upper bound for the number of critical points of the Green's function: 

Theorem 1.1. Suppose that the surface M is of finite type, that is, its fundamental 
group has finite rank. Then the number of critical points of any Green's function 
G on M is not larger than twice the genus of A4 , v, plus the number of ends, \, 
minus 1: 

^ critical points ^ 2;^ + A — 1 . 
// this upper bound is attained then G is a Morse function. 

In fact, the analysis of the set J- does not only yield this topological upper bound, 
but a more refined bound that exploits the conformal structure of the surface (see 



Theorem 5.1 ). This is particularly interesting because, as we shall see, it establishes 



some links between the conformal geometry of the surface and the portrait of the 
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gradient field VgG. It should be stressed that an analogous result does not hold 
for higher-dimensional Riemannian manifolds, as shown in |S]. 

A different proof of the estimates for the number of critical points of G, relying 
on methods from elliptic PDEs, was recently given by the authors in [^. However, 
the dynamical approach taken in the present paper provides a very satisfactory 
picture of the invariant sets connecting the different critical points of G and the 
dynamics of the field VgG, which cannot be obtained with the PDE techniques 
used in 6^. 

The article is organized as follows. In Section [2] we will present some basic 
facts regarding Green's functions on surfaces, including their obtention through an 
exhaustion procedure, their behavior at infinity and their connection with the con- 
formal structure of the surface. In Section [3] we describe the dynamics of the field 
VgG in a neighborhood of the pole y or a critical point. In SectionHlwe introduce a 
convenient compactification of the surface and establish some key properties of the 
sets D and T introduced in the Heuristic Principle above (and of some compactifica- 
tions thereof) . The structure of the set T and its compactification is characterized 
in Section [51 which allows us to prove the upper bound for the number of critical 
points of the Green's function. To conclude, in SectionlHlwe discuss the connection 
between the dynamics of the field VgG and the conformal geometry of the under- 
lying surface and make some comments regarding surfaces of infinite topological 
type. 

2. Green's functions on surfaces 

In this section we will recall what a Green's function is and how to obtain them 
using an exhaustion procedure, placing a special emphasis on conditions ensuring 
that the Green's function is "well behaved" at infinity, which is a key ingredient 
in the analysis of the dynamical properties of the vector field ^ gG. We will also 
discuss how to exploit conformal isometries to classify the possible behavior of the 
Green's function at the ends of the surface. Throughout this paper, the surfaces 
will be of class G°° and of finite topological type (that is, with finitely generated 
fundamental group) unless stated otherwise. 

The reason why it is crucial to make assumptions on the behavior of the Green's 
function at infinity can be readily seen even in the simplest case: the Euclidean 
plane M^. Indeed, if we let h be any harmonic function, it is clear that any function 
of the form 

Q{x,y) = -—\og\x-y\ +h{x) + h{y) 



is symmetric and satisfies the Green's function equation (1.1). The standard way 
of deciding which of these Green's functions should be "admissible" is to demand 
that the Green's function be obtained as a limit of Dirichlet Green's functions 
associated with an exhaustion of the plane by compact subsets (more details on 
this point will be given below). The only Green's function arising from such an 
exhaustion procedure would be the usual one, 

(2.1) e(x,y) = --^log|x-y|, 

ZTT 

which is the relevant Green's function for all geometric or analytical considerations. 
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In a general Riemannian surface {M,g), we will assume that the Green's func- 
tion G we consider shares the following two properties with the above Green's 



function (2.1 1 of the Euclidean plane. The first assumption is a weak monotonicity 
property for the Green's function on circles. The second assumption says that, 
when the surface admits a minimal (positive) Green's function (as is the case of the 
hyperbolic plane, although not of the Euclidean one), we should always consider 
this Green's function, for it plays a very special role in analysis and geometry. Re- 
call that a positive Green's function is minimal when it is pointwise smaller than 
any other positive Green's function. Observe that, if the surface (M,g) does not 
admit a minimal Green's function, then the infimum of any Green's function G over 
the surface is — oo. 

Assumption 1 (Monotonicity). The Green's function G is nondecreasing in the 
sense that 

(2.2) sup G = max G for all r > , 

M\Bg[y.,r) £>Bg(y,r) 

where Bg{y,r) denotes the geodesic disk in M centered at the pole y of radius r. 

Assumption 2 (Minimality). G is the minimal Green's function whenever the 
Riemannian surface admits a positive Green's function. 

It should be stressed that on any Riemannian surface (M, g) there are Green's 
functions that satisfy these assumptions. Indeed, the way one shows there always 
exist Green's functions on any surface is by taking a suitable limit of the Dirichlet 
Green's functions of an exhaustion of the surface by bounded domains. Since the 
Green's functions one obtains in this fashion actually satisfy the above assumptions, 
throughout this paper we will assume that the Green's function satisfies Assump- 
tions^ and^ For the benefit of the reader, we will next review this construction of 
Green's functions through an exhaustion by compact sets, which was introduced in 
the context of general Riemannian manifolds by Li and Tam |11| . To simplify the 
exposition, we will consider the case of Green's functions G(x) with a fixed pole y, 
but actually the procedure automatically yields the symmetric function Q(x,y). 

Consider an exhaustion fti C fl2 C ■ ■ ■ of the surface M by bounded domains. 
We can assume without loss of generality that the pole y belongs to the first do- 
main f2i. The idea now is to impose Dirichlet boundary conditions on each bounded 
domain flj and consider the corresponding Green's function Gj : nj\{y} — > M, 
which satisfies the equation 

AgGj = —5y in ilj , Gj = on dflj . 

One would be tempted to define the Green's function G as the limit of Gj as j — >■ oo. 
However, this limit does not exist in general. What can be proved (see e.g. [TT]) 
is that, for any choice of the domains flj, one can take a sequence of nonnegative 
real numbers {aj)°°^-^ such that Gj — Oj converges uniformly on compact sets of 
M\{y} to a Green's function G with pole y. Generally, the Green's functions 
obtained through this procedure are non-unique, but they satisfy the monotonicity 



assumption (2.2) and, when the surface admits a positive Green's function, this 
construction always yields the minimal one. 

Incidentally, it is worth pointing out that Green's functions do not exist on 
closed surfaces, which is the reason why we just consider noncompact surfaces in 
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this paper. To see why, it is enough to suppose that there is a solution of Eq. ( 1.1 ) 
in a closed surface M and, with a slight abuse of notation, integrate both sides of 
the equation over the whole surface and integrate by parts, which would yield the 
contradiction 

0= f AgG =- f 5p = -l. 

Let us now pass to describe how one can utilize conformal isometrics to analyze 
the behavior of the Green's function at each end. Recall that two Riemannian 
surfaces {M,g) and (A^,5o) ^tre conformally isometric if there is a diffeomorphism 
i" : A^ — >■ M and a smooth positive function f on A4 such that ^*g — fgo. An 
important property of the Laplace equation on surfaces is its conformal invariance, 
that is, if G{x) satisfies the equation 

AgG = -5y 

on the surface M for some point y, and another surface AA is conformally isometric 
to M through a diffeomorphism $ : A^ -^- A/, then 

G{x) := G(a>(x)) 

is a Green's function of A^ with pole y := $^^(j/): 

(2.3) A,„G = -^,. 

Furthermore, the corresponding gradient fields are orbitally conjugated through the 
relation 

(2.4) V,„G = /<f>*(V,G). 

A key result in the conformal geometry of surfaces, which will be of great use in 
this paper, is the uniformization theorem: 

Theorem 2.1 (Uniformization). There is a compact surface S of genus v with a 
metric of constant curvature go, a certain number Ai ^ o/ isolated points pi G E 
and another number A2 ^ o/ closed disks D^ G T, with smooth boundary such that 
the Riemannian surface M is conformally isometric to {A4,go), with 

(2.5) M:^J:\(\J{p,}u\Jd, 

^i=i j=i 

As is customary, we will call the deleted points {pi, . . . ,PAi} the parabolic ends 
of the surface M, while the deleted disks {Di, . . . , Dx^} are its hyperbolic ends. 
Furthermore, we will say that the parabolic end Pi is a removable singularity if the 
function G can be extended so as to satisfy the equation 

\oG = 

in a neighborhood oi pi in S. It should be noticed that an end being parabolic or 
hyperbolic is a geometric property of the surface, related to its conformal structure. 
However, whether a parabolic end is removable or not is not a geometric issue, as 
it depends not only on the surface M but also on the Green's function we consider. 

In the following two propositions we will relate the existence of parabolic and 
hyperbolic ends with the behavior of the function G at each deleted point or disk: 
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Proposition 2.2. // all the ends of the surface {M,g) are parabolic, the surface 
does not admit any positive Green's functions. Moreover, at each end pi we have 
that either 

(2.6) lim G{x) = -oo 

or Pi is a removable singularity. There is at least one point pi where the condi- 



tion (2.6) is satisfied. 



Proof. When the number of hyperbohc ends A2 is 0, it follows from Eq. (2.3) that 
the function G satisfies the equation Ag^G = everywhere in E but at the pole 
y and the isolated points pi. Furthermore, by Assumption 111 G is upper bounded 
at each point pi. If it is also lower bounded, it is standard that pi is a removable 
singularity [8] and A^gG = in a neighborhood oi Pi. If G is not lower bounded, pi 
is an isolated singularity of G, and the fact that G is upper bounded readily implies 
that AggG = CiSp. in a neighborhood oipi for some non- negative constant q. Hence 

_ ^^ 

(2.7) Ag,G^-Sy + J2c^Sp^ 

i=l 

in the closed manifold S, so integrating this equation over E and using that 
/j^AgpG = we infer that ^^Ci = 1. Hence Ci is positive for some i and, in 
view of the asymptotic behavior of any Green's function at a pole, it stems that G 
tends to —00 at the corresponding point pi. D 

Proposition 2.3. // the surface (M, g) has at least one hyperbolic end, there is 
a minimal positive Green's function G. The corresponding function G tends to 
zero at the boundary of each disk Di and all the parabolic ends pi are removable 
singularities. 

Proof. When the number of hyperbolic ends is A2 ^ 1, it is standard (for example, 
due to the existence of nonconstant positive harmonic functions in M [llj) that the 
surface admits a positive Green's function. Therefore, Assumption [2] ensures that 
G is the minimal Green's function of the surface, which corresponds to the unique 
solution G of the boundary problem 

A2 
Ag„G=-(5,j inEVljDj, G = on 91?^ for aU j. 

In particular, if there are any parabolic ends, they are all removable singularities. 

D 

Because of these propositions, one can extend the function G and the gradi- 
ent field Vg„G to all the removable singularities of the surface. We will find it 
occasionally convenient to consider this extension, which we will not distinguish 
notationally (it will be clear from the context). To conclude this section, we will 
present two examples that illustrate the issue of uniqueness and non-uniqueness of 
Green's functions on surfaces. 

Example 2.4. Consider the plane M^ with its Euclidean metric. It has a parabolic 
end, so that it is conformally isometric to the round sphere (§^,30) minus a point 
p via a diffeomorphism <I> : S^\{p} — > M? . 
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The standard Green's function with pole y, 

(2.8) G{x):=-^log\x-y\, 

obviously satisfies Assumption [T] Moreover, it is the only Green's function satis- 
fying this assumption, even though the plane does not admit any positive Green's 
functions. To see this, take the function G corresponding to any Green's function 
satisfying Ass ump tion fll Since the assumption is satisfied, Eq. (2.7) in the proof 



of Proposition 2.2 shows that G must satisfy the equation 



Ag„G ^ -Sy + Sp 



in the whole §^. Hence G is uniquely determined, so that the Green's function 



must be given by (2.8) 



Example 2.5. Let us consider the flat cylinder M x S^ and natural coordinates 
{z,6). It is conformally equivalent to the round sphere (S^,go) rainus two points 
{^15^2}, so it does not admit a positive Green's function. 
A Green's function with pole at (zq, 0o) is 

Gi{z,9) := log |"cosh(z — zq) — cos(6' — ^q)] . 

47r '- ■' 

This Green's function satisfies Assumption [l] It is not the only Green's function 

on the flat cylinder with this property; e.g., one can check that 

G2{z, 0) := -|^ log [ie^^ + ie^^" - e'+'° cos(0 - 9o)] 
is another instance. These Green's functions are connected by the identity 

G2{z,e)^G^{z,0)-'-±^. 

Notice that Gi tends to —00 at both ends (that is, as z — j' ±00) while G2 tends 
to — cx) as z — )- 00 but the end z -^ — cx) corresponds to a removable singularity 
ofGa. 

3. Local dynamical properties of Green's functions 

In this section we will carry out a local study of the dynamics of the gradient 
of the Green's function in a neighborhood of the pole or a critical point. Since the 
fields VgG and Vg^G are orbitally conjugated (cf. Eq. ( |2.4[ )), for convenience we 
will work with the latter gradient field instead. 

The first proposition we will prove here asserts that, when multiplied by a suit- 
able factor, the gradient vector field V^^G can be smoothly linearized at the point 
y, and that the corresponding normal form is a stable node. In particular, the 
trajectories approach the pole with a well-defined tangent. 

Proposition 3.1. There are C°° coordinates (xi,X2); defined in a neighborhood U 
of the point y in Ai and centered at y, and a smooth nonnegative function p : f7 — > M 
that only vanishes at the pole, such that the gradient field Vg^G can be written as 

d d 



PV,„G = -Xi — -X2^^^ 



in the domain U . 



2 - So ■ 
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Proof. Let us take local isothermal coordinates u = (ui, U2) centered at y, in which 
the metric reads as 

90 = f{u) (duj + duf) 
for a positive function /. Therefore, Eq. ( |2.3[ ) can be written in these coordinates as 

so G must be of the form 

G = -—\og\u\ + h{u), 
with |wp := u\ + ^2 and h a harmonic function: 

The gradient of G can then be expressed as 

with p{u) := 2iTf(u) |itp. The origin is then a hyperbolic zero of the vector field 
pVggG (which can be smoothly extended to the origin) and the corresponding 
eigenvalues are (—1,-1). Hence Siegel's theorem ensures that pWggG can be lin- 
earized via a diffeomorphism that is an analytic function of the coordinates (mi, U2) 
and the claim follows. D 

In the following proposition we will characterize the structure of the trajectories 
of the field ^ggG in a neighborhood of a critical point of G. In particular, we 
see that the trajectories approaching the critical point have a well-defined tangent. 
This proposition can be seen as a dynamical analog of Cheng's result on the critical 
points of eigenfunctions on surfaces [2]. Let us recall that the stable (resp. unstable) 
set of a zero z of a vector field is given by the points whose w-limit (resp. a-limit) 
is exactly the point z. 

Proposition 3.2. Let z be a zero of the vector field Vg^G (possibly a removable 
singularity pi) and let m ^ 2 be the degree of the lowest nonzero homogeneous term 
in the Taylor expansion of G — G{z) at z (which is always finite). Then z is an 
isolated zero and the intersection of a small neighborhood of z in Y, with either its 
stable or unstable set is homeomorphic to the set 

{CeC:r e[0,i)}, 

In particular, the index of the point z is 1 — m. 

Proof. Let us take isothermal coordinates u ~ {ui,U2) centered at z, in which the 
metric reads as 

50 = /(") {duf + dul) . 
Therefore, G is a harmonic function of u with respect to the fiat metric: 
.0 , d^G d^G „ 

It is therefore standard that G is an analytic function of {ui,U2). 
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Let hm be the first nonzero homogenous polynomial that appears in the Taylor 
expansion of G in these coordinates at 0, so that (with a slight abuse of notation) 

(3.2a) G{u) - G(0) = h.,^{u) + 0{\ur+^) , 

(3.2b) f{u) Wg„G{u) = Wh,n{u) + 0{\ur) , 

Here we are denoting by V the flat space gradient with respect to the coordinates u. 



By Eq. (3.1), the homogeneous polynomial hm is a harmonic function (with 
respect to the flat space Laplacian in the coordinates u) , so it must be of the form 

/i™ = CRe[e'"(ui+iu2)'"] 

for some real constants C, a. In particular, is an isolated critical point of hm, 
which readily implies that z is an isolated zero of V^qG. 

Let us now consider polar coordinates (r, 6*) € (0,e) x §^ defined by (mi,M2) — 
[r cos6,r sin6'). In these coordinates one has 

hm{r,0)^Cr"'cos{me + a). 

There is no loss of generality in setting a = 0. We define the polar blow up of the 
gradient Vg^G at z using polar coordinates as the vector field 



where we have used Eq. (3.2b). The blown- up trajectories are then given by 

(3.3a) f — r cosmO + O(r^) , 

(3.3b) ^ = -sin m6l + 0(r). 

The blown-up critical points are thus (0, 9k), with 6k '■= kn/m and fc = 1, . . . , 2m. 
The Jacobian matrix of X at {0,0k) is 

(3.4) DX{0,ek)=(^ ^~Q^' (_^*j,+i 

so these critical points are hyperbolic saddles. By blowing down, we immediately 
find that a deleted neighborhood of consists exactly of 2m hyperbolic sectors of 
the vector field X. 

Since the field X is proportional to the gradient field Vg^G through a factor that 
does not vanish but at z, this shows that the intersection with a small neighborhood 
of z with the its stable or unstable set is homeomorphic to 

{CeC:r e[0,i)}, 

as claimed. Besides, the well known Bendixson formula for the index of a planar 
vector field asserts that the index of z is 

, number of hyperbolic sectors 
md(z) = 1 = 1 — m , 

as claimed. D 

It is worth mentioning that the dynamics of the gradient of a harmonic function 
in a neighborhood of a critical point in dimension higher than 2 is much more 
involved, as is discussed in [S]. 
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4. The basin of attraction of the pole 

In this section we will introduce the concept of basin of attraction associated 
with Vgo G, as well as some convenient compactifications thereof. We shall see how 
this object and its boundary define a decomposition of the surface as the union of 
a disk and a 1-skeleton, as outlined in the Heuristic Principle in the Introduction. 
We shall also establish some properties of these sets. 

The basin of attraction of the pole y is a key object in this paper, and can be 
thought of as the set of points of the surface M that approach y when flowed along 
the trajectories of the field VgG. However, in view of the characterization of M 



in terms of a compact surface S (the Uniformization Theorem 2.1), it is slightly 
more convenient to study the basin of attraction directly in this compact surface, 
so instead we will use the following 

Definition 4.1. The basin of attraction is the set of points D in Ai whose w-limit 
along the trajectories of Vg^G is y: 

(4.1) D -.^ {x e M : uj{x) ^ y} . 

Of course, by the relationship between V g^^G and ^ gG^ the diffeomorphism $ : 
A^ -^ M maps the basin D onto the set of points in M whose w-limit along the 
integral curves of VgG is the pole y. It is easy to prove that D is diffeomorphic to 
a disk. 

It is clear that both D and its complement in AA are invariant sets under the 
local flow of V g^G. The complement, 

F:=M\D, 

will be a crucial object in the rest of the paper. 

In the following proposition we shall prove a general result about sets of M. 
that are invariant under the local flow of Vg^, G from which it stems an important 
property of J-" as a corollary: that M\D has empty interior and thus J^ coincides 
with the boundary of the basin D va M.. Notice that, since the basin of attraction 
is contractible, J- cannot be empty unless the surface M is diffeomorphic to M? . 

Proposition 4.2. Let S G M be an invariant set under the flow of ^ g^G that 
does not contain the point y and is relatively closed in M.. Then the interior of S 
is empty. 

Proof. In this proof, we will assume that we have enlarged the set A^ and extended 
the function G in the obvious way so that the removable singularities are points 
contained in Al and G is also defined at these points. Let U denote a connected 



component of the interior of S. By Propositions 2.2 and 2.3 the disks Di or 
the deleted points pi that are not removable singularities behave as local minima 
of the function G. Since ^g^G = both in U and in a neighborhood of any 
removable singularity pi, the maximum principle for harmonic functions ensures 
that the maximum of G must be attained at a point z of dU (possibly a removable 
singularity) . 

Since the boundary of U is inva rian t, z must be a critical point of G, which is 



necessarily isolated by Proposition 3.2 As G is smooth in a neighborhood of z, z 
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is an isolated maximuiii of G in the closure U and G is increasing along the local 
flow V't of Vgp G, it follows that for any e > there exists some (5 > such that 

^lJt{Bg,{z,5)f^U)(lBg,{z,e) 

for alH > and 

Hence there exists a region Bg^{z,S) (lU oi nonzero measure whose w-limit is z. 
This contradicts the fact that, since ^ggG = in a neighborhood of z, the local 
flow of Wgg G is area-preserving, so the set U must be empty. D 

Hence we immediately obtain the desired statement about T: 

Corollary 4.3. The complement of D in Ai coincides with the boundary of D in 
A4, that is, 

T = M\D = dD . 

Our goal now is to derive further properties of the set T . For technical reasons, 
it is more convenient to consider the flow of an auxiliary vector field X defined in 
the whole compact surface S rather than that of V g^G, which is only defined on 
M.. For this, let us take a point qi belonging to the interior of each disk Di. Since 
the disk Di retracts into qi, one can take a diffeomorphism 

■^ : M-^ Y\{pi,...,px^,qi,...,q\^} 

which is equal to the identity outside a small neighborhood of the closed disks Di 
(in particular, at y). 

Let us relabel the parabolic ends if necessary so that {pi}i^i are the non- 
removable singularities, with ^ A'j^ ^ Ai. The auxiliary field is then defined 
as 

X:=F^,{Vg„G), 
where F : S — > M is a smooth nonnegativc function that only vanishes at {y} U V , 
with 

(4.2) "P := {pi,...,pA;,gi,---,9A2} 

and is chosen so that X can be smoothly extended to the whole surface E. It is 
standard that such a factor always exists; notice, moreover, that we do not need to 
impose any additional conditions on this factor to ensure that X is also defined at 
the removable singularities (we only need to consider the obvious extension of G). 
Throughout this paper, we will denote the flow of X hy (pt- 

It is clear that the function Go$^^ can then be extended to a continuous function 
G : S — > [— oo, +oo] by setting 

G(y):=+oo, G{pi) := -oo , G((Z,) := 

for 1 ^ i ^ A']^ and 1 ^ j ^ A2. Hence it stems that the field X is gradient-like: 
the Lie derivative CxG is strictly positive but at the zeros of the field X, which 
are the only points at which the function G can fail to be smooth. Moreover, 
the zeros of X are exactly the images under the diffeomorphism ^ of the zeros of 
Vg(jG (including those that correspond to removable singularities) and {y} U V. 
Notice that the points in V are precisely the global isolated minima of G: indeed. 
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by Proposition 2.3 if A2 ^ 1, then G is positive and X[ = 0. Besides, from this 



proposition it stems that the cardinahty of V is 
(4.3) A' := 




The reason why we consider the field X is that, instead of directly analyzing the 
sets D and T associated with the field Vg^G, it is easier to consider the analogous 
dynamical objects for the field X. That is, we will consider the set D of the points 
in E whose w-limit along the flow of X is y, which we will still call the basin of 
attraction of the field X. Its associated basin boundary is then defined as 

^:=S\i5. 

Just as in the case of the set D, it is standard that D is diffeomorphic to a disk. 



Remark 4.4. Proposition |4.2| and its proof are also valid, mutatis mutandis, for the 
field X. That is, if 5 C S is a closed invariant set under the flow of X that does 
not contain the point y, it has empty interior. In particular, 

^ = j:\d = dD . 

The following proposition completely characterizes the set T in terms of the zeros 
of the field X and their stable components. To state this result, let us introduce 
the notation 

(4.4) C:=^[{zeMU {pv,+i, • • ■ ,PaJ : Vg„G(z) = O}] 

for the image under the diffeomorphism '5 of the critical points of G, including 
those that may correspond to a removable singularity. 

Proposition 4.5. J^ is the union of the zeros of the field X other than y and the 
stable sets of the zeros in the set C: 

zGC 

Furthermore, J- is connected. 



Proof. Remark 4.4 readily implies that T is connected and has empty interior. Let 
us now set 

W -.^VulJ W'{z). 
zee 
By definition, it is clear that the w-limit of any point x d W cannot be the point 
y, so T^ C ^\D. 

Let us now prove the converse implication: S\_D C W. For this, we shall show 
that the w-limit of any x G I]\Z? must be a zero of the field X different from y. We 
can obviously suppose that x is not a zero of X. Since X is gradient-like, (f>tx must 
tend to a zero z oi X as t —> 00 [15 . Besides, z must belong to the set C because 
the points in V are minima of G and the Lie derivative 

CxG{cj)tx) = ^d{cj)tx) 
at 

is positive if x is not a zero of X. Hence we infer that !F = W. D 
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It is clear that a good description of T immediately yields a complete charac- 
terization of the original set T . In particular, Proposition |4.5| easily implies that, 
roughly speaking, T consists of the set T and some additional segments that con- 
nect a point in T with a parabolic end of the surface. From this set, of course, we 
still have to remove the points in J- corresponding to an end. This is the content 
of the following 

Corollary 4.6. The set T is the image under the diffeomorphism ^~^ of the set 
-?U y {(j)_tPi-t>Q}\\{pi,...,Px^,qi,...,q\^} . 

The closure of J- in T, is connected. 

Proof. It is clear that the image under vj/^i of the set J- minus the ends 

{pi,...,px^,qi,...,qx^} 

must be contained in J-'. Likewise, the image under ^~^ of any point x in D will 
be in the basin of attraction D unless at some positive time t the trajectory (j>tx 
passes through an end, which is necessarily a removable singularity since the other 
ends are all zeros of the field X. That is, for some i > one would have 

X = (j)-tPi for some A'j^ -I- 1 ^ i ^ Ai . 

This proves the formula in the statement. (Of course, among the removable singu- 
larities it would be enough to consider those that are not a zero of the field X.) 

Therefore, the closure of J^ in S is either empty (then F consists of a single 
point) or is diffeomorphic to the union of J- and a finite number of segments whose 
endpoints are a removable singularity pi as above and a zero of the field X (which 
obviously belongs to F). Therefore, the closure of F is connected. D 

5. Structure of the basin boundary and bounds for the critical 

POINTS 

In this section we will provide a full characterization of the basin boundary, 
thereby obtaining an upper bound for the number of critical points of the Green's 
function. The characterization of the basin boundary lay bare a strong connection 
between the dynamics of the field X (or, equivalently, V gG) and the topology and 
conformal structure of the surface. 

The following theorem provides an upper bound for the number of critical points 
of G (including those corresponding to removable singularities) in terms of the 
conformal properties of the surface, which appear through the number A' introduced 



in Eq. (4.3). A straightforward consequence of this result is the purely topological 
bound presented in Theorem |1.1[ which differs from the present statement in that 
here we are using additional information on the conformal structure of the surface 
S to sharpen the upper bound: 

Theorem 5.1. The number of zeros of the field V^qG in 

7Wu{pa;+i,---,PAi}, 
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that is, the cardinal of the set C , is not larger than 2i/ + A' — 1, and if this upper 
bound is attained then G is Morse. 

Before presenting the proof of this result, it is illustrative to sketch the argument 
in the easiest case: when M is diffeomorphic to M? . In this case, the proof is similar 
to that of a result on the absence of critical points in some boundary value problems 
in the exterior of a bounded domain in M" that we proved in 5]. Of course, this 
particular case is elementary and could be easily treated using the uniformization 
theorem, but it serves to illustrate the basic dynamical ideas underlying the proof 
of the general situation. 

Sketch of the proof when M is diffeomorphic to M? . Let us analyze what happens 
in M , which is diffeomorphic to the sphere minus a point p. Suppose that we have 



a zero z. By Proposition 3.2 its stable set (with the point z deleted) consists of at 



least two curves. The a-limit of each of these curves cannot be y by Proposition 3.1 
so as the vector field X is gradient-like either the curve approaches the end p or its 
a-limit is another zero zi. We can now apply the same argument for the zero zi, 
and if necessary to successive zeros 2:2, Z3, . . . Notice that Zj 7^ zj. for j 7^ k, since 
otherwise we could have an invariant set (defined by the union of curves in the 
stable sets of zeros of the field) that is a Jordan curve not containing y. Then this 
invariant set would separate the plane in two disjoint invariant sets with nonempty 



interior, contradicting Proposition 4.2 



Since the zeros are isolated in A4 by Proposition 3.2 we can eventually take a 
union of curves in the stable sets of zeros Zj of the field (possibly infinitely many, but 
only accumulating at p) whose closure in the sphere is a Jordan curve that contains 
the point p but not y. Again, this curve encloses an invariant set with nonempty 
interior that does not contain the point y, in contradiction with Proposition |4.2[ 
Hence there cannot be any zeros of '^g„G and the theorem follows. D 

In the general case, the proof is more involved and relies on a careful analysis of 
the saddle connections between zeros of the field. In the demonstration we need two 
lemmas that are presented right after the proof and make use of the same notation. 

The proof is divided in two parts. First we show that the number of zeros of X is 
finite. Notice this is not trivial, since the zeros of X could accumulate at the set V, 
which is associated with ends of the surface M , wit hout contradicting the fact that 



critical points of G are isolated by Proposition 3.2 However, we show that if there 
were an infinite number of critical points, there would be infinitely many closed 
invariant curves in S defining independent homology classes, which is forbidden by 
the fact that the fundamental group of S has finite rank. Roughly speaking, these 
closed invariant curves are constructed by successive continuation of the stable sets 
of some zeros of X. The second part of the proof consists in estimating the number 
of zeros of X using Hopf 's index theorem and the characterization of the dynamics 
of X in a neighborhood of each zero. 



Proof of Theorem \5.1\ Let us start by recalling that the set V, introduced in Eq. (4.2 1, 
consists of precisely A' points, which are the minima of the function G. As X is 
gradient-like, the a- and w-limit sets of any trajectory of this field are necessarily 
a zero of X. Moreover, the zeros of this field are obviously given by the set 

Z:^{y}\JC\JV, 
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where we record here that the set C, defined in (4.4), corresponds to the critical 



points of G (possibly including removable singularities) under the diffe omorphism 



^ . The fact that the critical points of G are isolated (by Proposition 3.2 ) guarantees 
that Z does not accumulate but possibly at V . 

Let 7 be a trajectory of the field X. This trajectory will be called constant 
if it consists of a single point. We will use the notation a(7) and w(7) for the 
a- and w-limit sets of 7. Let us introduce a partial order on the set of zeros Z 
as follows. Given two points x,x' € Z, we shall write x ^^ x' if, for any open 
neighborhoods U 3 x and y 9 x' in E, there exist integers p ^ 0, g ^ and 
nonconstant trajectories (of the field X) 7j,, . . . , 7g such that 

(i) uj{jp) e U, a(7,) e V, 
(n) a{-fj) = w(7j+i) for p s$ j < g 



1. 



Finiteness. We claim that the set of zeros Z is finite. In order to prove this, let us 
assume the contrary. By Lemma |5.3| below, for each point a; e C there exists some 
point p Cz V such that x )~ p. Since V is finite, there exists some pi € V such that 
one can choose a sequence {xk)^^i of distinct points in C with Xk >- P- For each 



point Xk, Lemma 5.2 below yields a continuous path Ffe.i : [0, 1] — > S whose image 
is invariant under the field X and satisfies rfc^i(O) = Xk and rfe,i(l) = pi. 

As a straightforward consequence of Proposition |3.2[ the invariant set 

w%xk)\rkA[o^]) 

is nonempty. If we let Xk be the a-limit of a trajectory contained in this set, 
we obviously have x^ >- x^- By Lemma |5.3[ either Xk & V ot Xk > p for some 
p G V. Since V is finite, by Lemina |5.2| and possibly upon restricting ourselves to a 
subsequence that we still denote by (xfc)^]^, we obtain a family of continuous paths 
rfc,2 : [0, 1] -T> S whose image is invariant under X and such that rfc_2(0) = Xk and 
rfc,2(l) = P2 for some fixed P2 GV (possibly the same as pi). 
By construction, for each positive integer k the connected set 

rfe,i([o, 1]) u rfe^2([o, 1]) u rfc+i,i([o, i]) u r,+i,2([o, i]) 

consists of a continuous curve that connects the points x^ and Xk+i passing through 
Pi and another continuous curve that connects the same pair of points Xk,Xk+i 
passing through p2- It is then evident that this set, which can have self-intersections, 
contains an invariant loop (continuous closed curve) A^. One can obviously ensure 
that Afe ^ Ak' for k ^ k' (these loops can intersect, though), and that the point y 
does not belong to any A^. 

For any integer j, the union of invariant loops 



UA. 



fc=i 

cannot disconnect S, since a connected component of S\lJ^^-^Afe that does not 
contain y would be an invariant set with nonempty interior, contradicting Propo- 
sition 4.2 Therefore, it is standard that the homology classes [A^] € Hi{I];Z) 
defined by the cycle A^ must be independent for all fc = 1, 2, . . . This is impossible 
in a finitely generated surface, so we infer that the set C is finite. 
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Upper bound. Let us now pass to bound the cardinality of the set C. Suppose that 
the number A2 of removed disks is at least one, so that 

V = {qi,...,qx^}. 

A first observation is that the points qi are isolated zeros of the field X, by the 
finiteness of Z, and are local repellers because they correspond to minima of the 
function G. Therefore, the index of X at these points is 

ind(gi) = 1 . 

Similarly, the point y is an isolated zero which is a local attractor, so it has index 
1. If we now apply Hopf's index theorem to the vector field X in E, we get that 
the sum of the indices of the zeros of X equals the Euler characteristic of S: 

A2 
(5.1) ind(y) + J2 ind(^) + £ ind((Z,) = x(S) = 2 - 2iy . 

zee i=i 



Since the index of each point z e C is smaller than or equal to —1 by Proposition 3.2 
plugging the values of the indices of y and qi we find 

(5.2) #C sC -^ind(z) = 2!y + A2-l. 

zee 
The equality is not satisfied but perhaps when ind(z) = — 1 for all z G C, that is, 



when G is Morse (by Proposition 3.2 1. This proves the theorem when A2 ^ 1. 
Consider now the case where A2 = 0, so that 

V^{pi,...,px'^}. 
Arguing as before one easily finds that the index of X at each point pi is 

ind(pj) = 1 . 
We now apply Hopf's index theorem to the vector field X in S to find 

(5.3) ind(y) + ^ ind(z) + ^ ind(pO = x(S) = 2 - 2^. , 

zee i=i 

so that the same argument as above yields 

(5.4) #C < - ^ ind(z) = 2j/ + A'l - 1 . 

zev 

Again, the inequality being saturated at most when ind(z) = —1 for all z G C (that 
is, when G is Morse). The theorem then follows. D 

Lemma 5.2. Let x,x' Cz Z such that x )~ x' . Then there exists a (not necessarily 
unique) injective continuous path F : [0, 1] — > S such that: 

(i) r(0) ^x andr{l) ^ x' . 
(ii) G oT is strictly decreasing. 
(iii) The curve r([0, 1]) is invariant under the field X . 

Proof. By the definition of the partial order >- and Zorn's lemma, there exists a 
countable sequence {7j}'^^p {—p,q S N U {00}) of nonconstant trajectories of the 



field X satisfying Conditions (i) and (ii) in the proof of Theorem 5.1 and such that 



lirn a;(7j) = x , lira q;(7j) = x' . 
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Let us consider any continuous parametrization 



r:[o,i]^U^^Wc^ 



3=P 

mapping to a; and 1 to x' . Since the Lie derivative CxG is positive in T\Z, G 
is increasing along nonconstant trajectories, which inipUes that G o F is strictly 
decreasing (notice that the definition of F accounts for the fact that this function 
is decreasing instead of increasing). Therefore, F is injective. Moreover, the curve 
F([0, 1]) is clearly invariant because it is the union of trajectories. D 

Lemma 5.3. For any z €z C U {y} there exists some p G V such that z )~ p. 
Moreover, there are no x £ Z such that x )~ x or p )~ x for some p € V (i.e., the 
elements in V are minimal with respect to the partial order). 



Proof. Let us take an element z e C. By Proposition |3.2[ there exists a neighbor- 
hood U of z such that {^{z) n U)\{z} has at least two components Ci,C2 and 
each Ci is a piece of a trajectory of X. Since X is gradient-like, the a-limit set of 
Ci is another zero zi E CUV. 

If zi G V, the statement follows. Otherwise, we can repeat the previous argument 
replacing z hy Zi. As X is gradient-like, proceeding this way we obtain a sequence 
of distinct points zi, Z2, . . . in Z. If there is some point Zk & V, we are done, so we 
can assume that there is a sequence (zfc)^^! of distinct points o f Z with Zk >- Zk+i- 
Since ^~^(C) consists of isolated points in Ai by Proposition 3.2 it follows that 
distgg{zi;,V) tends to zero as fc — > oo. Hence there must exist some p (£ V such 
that a subsequence {zk')k'eJ tends to p. As above, it then follows that the initial 
sequence {zk)'^i also tends to p because of the fact that p is an isolated minimum 
of G and G is increasing along trajectories. 

Notice that the proof also applies when we start with the point y instead of a 
point z (£ C. Finally, note that obviously x'lf x for any a: S 2^ by Lemma |5.2[ and 
that p )f X because each p G T' is an isolated minimum of G and G is increasing 
along the flow of the field X. D 

In the following theorem we show that the basin boundary J- encodes the topol- 
ogy of the surface E. In particular, the flow of the field X defines in a natural way 
a decomposition of S into a disk D and its 1-skeleton T . This decomposition is 
similar to the one arising when one considers the cut locus of a point in E (with re- 



spect to some metric go) but it is generally different. Together with Proposition 4.5 
and Corollary |4.6[ this provides a rigorous reformulation of the Heuristic Principle 
stated in the Introduction. 

Theorem 5.4. The set J- is a connected graph with the same homology as the 
surface T,: 



Proof. We showed in Proposition 4.5 that T consists of the zeros of X other than y 
(that is, CUT') and their stable components, which are continuous curves with 
endpoints belonging to C U 7^. Since C is finite by Theorem [O] it then follows that 
J^ is a connected graph. 
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Let B he a, small disk in S centered at y and let us consider its image under the 
time-t flow of X, 4>t{B). It is apparent that S\{y} deform retracts onto the set 

S\(/'t(S), 
for any i ^ 0. Moreover, 

^=n(sw*(s)) 

by the definition of the set J- . This ensures that J^ is a strong deformation re- 
tract of I]\{y}, so it is well known that, J- being a connected graph, we have the 
isomorphism of homology groups 

iJi(^;Z) = iJi(I]\M;Z). 

Since 

i/i(E\{y};Z) = Hii^\B;Z) = i/i(I];Z) = Z^'' 

by a standard argument using the Mayer-Vietoris sequence, the theorem follows. 

D 



Remark 5.5. Because of the connection between T and T, Theorem 5.4 also gives 
very detailed information about the structure of the set J^. In particular, J^ is a 
graph but is necessarily noncompact and possibly disconnected. Moreover, the rank 
of Hi{F]'L) is at most 2v but can be strictly smaller than this number, as some 
of the cycles that appear in J- can be killed after removing the points of E that 
correspond to the ends of the noncompact surface. 

To conclude, let us present an illustrative example in which the different sets 
that we have been discussing in Sections HI and [5] can be computed explicitly: 

Example 5.6. Let M be the torus minus the point p := (0, tt), written in terms 
of the standard 27r-periodic coordinates on the torus. We choose a complete, con- 
formally flat metric g on M, and flx the position of the pole at y = (0,0). It is 
then standard that Li and Tam's procedure [TT] gives rise to a Green's function G 
invariant under the isometric transformations of the flat torus that fix the point y. 

It is then clear that the curves 

{0} X S^ , {tt} X §1 , §1 X {0} , §1 X {tt} 

are then invariant under the local flow of VgG. An easy argument using the sym- 
metries then show that the points (7r,0) and (7r,7r) are zeros of the field VgG, 
which must be nondegenerate (hence hyperbolic saddles). These are the only zeros 
because the upper bound in Theorem |5.1| is attained. 

The set T consists of the two circles {tt} x S^ and S^ x {tt}, which respectively 
correspond to (two) saddle connections and to the trajectories of X that connect 
a saddle with the point p. Since J^ contains two independent cycles, Hi{T]'L) is 
isomorphic to the first homology group of the torus. The set F is then given by 

F = F\{p), 

so is consists of a closed curve and an open curve. In particular, Hi{J-]'L) = Z. 
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6. Applications and remarks 

In this last section we will make some remarks about the connections between the 
critical points of the Green's function and the conformal properties of the surface. 
We will end with some comments about Green's functions on surfaces of infinite 
topological type. 



Conformal structure and dynamics. A consequence of Theorem |5.1| is that by 
analyzing the dynamics of the gradient field V gG (or, equivalently, of the field X) we 
can sometimes extract information about the conformal structure of the underlying 
surface. For example, a surfaces must satisfy very stringent geometric conditions 
in order to admit a Green's function without any critical points, as we show in the 
following 

Proposition 6.1. Let us suppose that the surface M has at least two ends. If 
all its ends are either hyperbolic or non-removable singularities, then any Green's 
function on M has at least one critical point. 



Proof. Eqs. (5.2) and (5.4 1 mean that, with the same notation as in the proof of 
Theorem 15.11 



(6.1) -^ind(2) = 2t/-l + A'. 

zee 
Let us call 

Ci ■.^C\{px'^+i,...,px,} 
and C2 '.= C\Ci. Any of these sets can be empty, and it is clear that the cardinality 
of Ci equals the number of critical points of the Green's function G in M . The 
cardinality of C2 is at most Ai — A'j^. 

As A' — X[ + A2 , Eq. ( |6.1| can be rewritten as 

(6.2) - Y^ ind(z) = 21/ - 1 + a; + A2 + ^ ind(z) . 

zsCi zeC2 

If all the ends are hyperbolic, Ai — X[ = 0, and if all ends are non-removable 
singularities, then Ai = X'l and A2 = 0. Therefore, in both cases C2 is empty, so 



the RHS of (6.2 1 is nonzero and Ci cannot be empty. Since hyperbolic ends and 



non-removable singularities cannot coexist, the proposition is proved. D 

Surfaces of infinite topology. It is worth emphasizing that the case of surfaces 
whose fundamental group is not finitely generated is totally different from the case 
of surfaces of finite type. We shall next provide examples of surfaces that are not 
finitely generated both with an infinite number of critical points and without any 
critical points. 

Example 6.2. Let us consider the case where M is a torus of infinite genus with 
two refiection symmetries (see Figure IT]). For this, we can regard M as a surface 
embedded in M'^ and invariant under the reflections 

Ili{xi,X2,X3) -.^ {~xi,X2,X3) and 113(2:1, 0:2, ^3) := (xi, a;2, -xs) . 

We will endow M with the metric g induced by the Euclidean metric in M.^ and 
consider a Green's function G on M (satisfying Assumptions 1 and 2 in Section [2]) 
with a pole at a point y invariant under the reflections 111,113. Li and Tam's 
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Figure 1. A torus of infinite genus with two reflection symmetries. 



procedure can be used to obtain a Green's function with the same symmetries 
(that is 

G = G o Hj 
forj-1,3). 

It is straightforward that the intersection of the surface M with the plane {xj = 
0} is invariant under the local flow of ^ gG, for j = 1,3. A simple argument then 
shows that the intersection of M with the X2-axis must be invariant too. Since it 
consists of isolated points, these must therefore be either critical points of G or the 
pole y. These critical points are obviously infinite in number. 

Example 6.3. Let us consider the unit disk 

©2 :={(a;i,X2) : xl + xj < 1} 

with its hyperbolic metric gi. We let M be the unit disk with infinitely points 
removed as 



M = 



\((o,o)uU(o. 



1 



n=2 



The hyperbolic metric gi is not complete on M, but it is well known [THj that there 
is a conformally equivalent metric g = xQi such that (M, g) is complete. 

Consider the minimal Green's function G with a pole at a point y e M. By 
the conformal invariance of the Laplacian, this Green's function is precisely the 
minimal Green's function of the disk D'^ with the hyperbolic metric gi (or rather 
its restriction to M). Since the gradient of the minimal Green's function of the 

(Af , g) is an example of a 



5.1 



hyperbolic disk {D ,gi) does not vanish by Theorem 
surface that is not finitely generated whose minimal Green's function does not have 
any critical points. Incidentally, notice that, although the proof of Theorem 5.1 
does not apply to surfaces of infinite topological type, zero is exactly the upper 
bound one gets by recklessly applying the formula 2i^ — 1 + A2 in this case (notice 
that for M we have t^ = 0, Ai = 00, A'j^ = 0, A2 = 1). 



Acknowledgments 

This work is supported in part by the Spanish MINECO under grants FIS2011- 
22566 (A.E.), MTM2010-21186-C02-01 (D.P.-S.) and the IGMAT Severo Ochoa 
grant SEV-2011-0087 (A.E. and D.P.S.), and by Banco Santander-UCM under 
grant GR35/10-A-910556 (A.E.). The authors acknowledge the Spanish Ministry 
of Science for financial support through the Ramon y Cajal program. 



References 

1. M. Brelot, G. Choquet, Espaces et lignes de Green. Ann. Inst. Fourier 3 (1951) 199-263. 

2. S.Y. Cheng, Eigenfunctions and nodal sets, Comment. Math. Helv. 51 (1976) 43—55. 



GRADIENT DYNAMICAL SYSTEMS ON SURFACES AND GREEN'S FUNCTIONS 21 

3. S.Y. Cheng, S.T. Yau, Differential equations on Riemannian manifolds and their geometric 
applications, Comm. Pure Appl. Math. 28 (1975) 333-354. 

4. T.H. Colding, W.P. Minicozzi II, Large scale behavior of kernels of Schrodinger operators, 
Amer. J. Math. 119 (1997) 1355-1398. 

5. A. Enciso, D. Peralta-Salas, Critical points and level sets in exterior boundary problems, 
Indiana Univ. Math. J. 58 (2009) 1947-1968. 

6. A. Enciso, D. Peralta-Salas, Critical points of Green's functions on complete manifolds, J. 
Differential Geom. 92 (2012) 1-29. 

7. A. Gabrielov, D. Novikov, B. Shapiro, Mystery of point charges. Proc. Lond. Math. Soc. 95 
(2007) 443-472. 

8. D. Gilbarg, J. Serrin, On isolated singularities of solutions of second order elliptic differential 
equations, J. Anal. Math 4 (1955/56) 309-340. 

9. P. Goldstein, Gradient flow of a harmonic function in R"^, J. Differential Equations 247 (2009) 
2517-2557. 

10. A. A. Grigoryan, L. Saloff-Coste, Dirichlct heat kernel in the exterior of a compact set, Comm. 
Pure Appl. Math. 55 (2002) 93-133. 

11. P. Li, L.F. Tam, Symmetric Green's functions on complete manifolds, Amcr. J. Math. 109 
(1987) 1129-1154. 

12. P. Li, L.F. Tam, J. Wang, Sharp bounds for the Green's function and the heat kernel. Math. 
Res. Lett. 4 (1997) 589-602. 

13. P. Li, S.T. Yau, On the parabolic kernel of the Schrodinger operator. Acta Math. 156 (1986) 
153-201. 

14. C.S. Lin,C.L. Wang, Elliptic functions. Green functions and the mean field equations on tori, 
Ann. of Math. 172 (2010) 911-954. 

15. W. de Melo, J. Palis, Geometric theory of dynamical systems, Springer, Newr York, 1982. 

16. M. Morse, Equilibrium points of harmonic potentials, J. Anal. Math. 23 (1970) 281-296. 

17. M. Morse and S.S. Cairns, Critical point theory in global analysis and differential topology, 
Academic Press, New York, 1969. 

18. K. Nomizu, H. Ozcki, The existence of complete Riemannian metrics, Proc. Amer. Math. Soc. 
12 (1961) 889-891. 

19. R. Shelton, Critical points of harmonic functions on domains in R^, Trans. Amer. Math. Soc. 
261 (1980) 137-158. 

20. J.L. Walsh, The location of critical points of analytic and harmonic functions, AMS, New 
York, 1950. 

Instituto de Ciencias Matematicas, Conse.io Superior de Investigaciones Cienti'ficas, 
28049 Madrid, Spain 

E-mail address: aenclsoOicmat .es 

Instituto de Ciencias Matematicas, Consejo Superior de Investigaciones Cienti'ficas, 
28049 Madrid, Spain 

E-mail address: dperaltaOicmat .es 



